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We study the recently reported new characteristic gapless charge ordered state in a spinless 
fermion system on a triangular lattice under strong inter-site Coulomb interactions. In this state 
the charges are spontaneously divided into solid and liquid component, and the former solid part 
aligns in a Wigner crystal manner while the latter moves among them like a pinball. We show that 
such charge ordered liquid is stable over a wide range of filling, 1/3 < n < 2/3, and examine its 
filling dependent nature. 

PACS numbers: 



INTRODUCTION 

Ever since the concept of frustration-induced "spin liq- 
uids" is introduced to the triangular lattice systems yj, 
the theoretical search for a special state without any 
types of long range ordering has been proceeded in 
the the strongly correlated systems under geometrical 
frustration^, Q- Such liquids are now beingexplored 
experimentally in the materials like NiGa2S4[4| and k- 
ET 2 Cu 2 (CN) 3 |3. However, the competition among sev- 
eral orders towards their destruction often makes the sys- 
tem reconcile with another different type of ordering. 
Simplest examples are the spins on frustrated lattices 
forming dimersjg, plaquettesQ, etc. Most recent ones 
are reported on hard core bosons of a nematic ordering 
in frustrated square lattice @| as well as of supersolids in 
the triangular lattice 0]. We also proposed a frustration 
induced novel charge liquid in a triangular lattice denoted 
as a "pinball liquid" .10], which turned out to have very 
similar character with the above supersolids despite the 
difference in the statistics. This state has a coexistence 
of a Wigner crystal-type of solid and itinerate liquid type 
of charges so that a phase separation is expected with- 
out any domain structures. In this paper we examine 
the filling dependent nature of this liquid in detail in the 
spinless fermion system and confirm the strong coupling 
picture. 



MODEL 

The present spinless fermion system is described by the 
t-V model hamiltonian, 

H = ^ ( - tcjcj + h.c. + Vmn^j , (1) 

where cj denotes the annihilation operator of fermions 
and rij(= CjCj) is its number operator. The summation is 
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FIG. 1: (a) Schematic illustration of the representative clas- 
sical ground states, (a) two-fold stripes, (b) basic three sub- 
lattice structure and one of its contingent. 



given over the nearest neighbor(nn) pair sites on a trian- 
gular lattice with indices (ij). Together with the strong 
coupling argument, we perform the exact diagonalization 
on N = 4 x 6 cluster with periodic boundary conditions, 
which could fully describe both the three fold and two 
fold type of periodicity which reside in the present sys- 
tem. Since the size dependence of the energy is confirmed 
to be very small and also the correlation function shows 
a well converged character within a relati vely short range 
in comparison with the parallel works |9Ul Of. the present 
results are enough qualified to support the physical pic- 
ture we present. 



STRONG COUPLING THEORY 

In the strong coupling limit, t = 0, the system is iden- 
tical with the classical Ising S=l/2 spin system on a tri- 
angular lattice. The ground state is disordered due to 
macroscopic number of degenerate states 12] which are 
classified into two groups according to their periodicity 
as shown in Fig. ^ one portion is composed of two-fold 
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FIG. 2: Kinetic energy gain, E x = (Eq - E c )/N, in the 
formation of a pinball liquid at N — 24 with t — 1 as a 
function of (a) l/V and (b)iVe, where N e =8,12 corresonds to 
n=l/3, 1/2, respectively. The fitting by the power of 1/1/ 
is given in (a), whose results are reflected in (b) as those of 
V = oo. The band energy gain of the free electrons on a 
corresponding honeycomb(iV = 16) and triangular (N — 24) 
lattices are shown together for comparison. 




FIG. 3: (a) Example of two basis that is allowed to mix by t. 
The ball can move from site i to j when the extra two neigh- 
bors of site i with big circles together have the same number of 
balls with those of the site j. (b) Representative basis of the 
three sublattice states at n = 0.5. Filled black circle and gray 
or red circles denote the "pin" and "ball", respectively, and 
the open circle the empty C-sites. The "balls" on C-site can 
move to neighboring B-site and vise versa when condition (a) 
is fulfilled. Arrows show the possible movement, (c) Example 
of the case when the "pin" (big black circle) on A-site move 
to its neighbor. When the black circle moves away from the 
configuration shown in the bottom, the two "balls" marked 
with big circle can no longer move. If one of these "balls" 
refill the original "pin", then other balls can move away. 

striped states which have staggered alignment of charges 
in any one of three directions but is disordered in other 
two directions. This state has an order V^/V-degeneracy. 
The others have a three sublattice structure represented 
by Fig. IHb). Here, as long as A and C sublattices are 
filled and empty, respectively, the configuration of the B- 
sites is chosen arbitrary. In addition, there are large num- 
ber of contingents which include the irregular structure 
on A and B sublattices as shown in Fig. dc) but still has 
a same binding energy. Including all of them, the three 
sublattices can afford any fillings over 1/3 < n < 2/3 



FIG. 4: (a) three-body correlation function, Ps(i) = (ni(l — 
n2)rii), along the second chain at several fillings between 
1/3 > tl / 1/2. N e = 8,12 correspond to 1/3 and 1/2- 
fillings, respectively, (b) Amplitude of P3 as a function of N e 
at sites i — 4, 5, 6 which belong to the A, B, C-sublattices, 
respectively. The difference between pins and balls become 
small near half-filling. 

without changing the binding energy, in sharp contrast 
to the strictly half-filled stripes. 

The introduction of t 7^ drives the classical disor- 
dered system to a new ordered phase which we call a pin- 
ball liquid^. For simplicity, we start from the n = 1/3 
case which has a well known unique ground state where 
only the A-sublattice is filled. Then the extra binding 
energy 2V is required to move a single charge from one 
of the A-sites, thus the "pins" on A-site form a solid. 
Next we add a single extra charge on one of the empty 
sites. This charge can move to another nearest empty 
site without losing the binding energy. Therefore once 
the filling is off n = 1/3, the system discontinuously be- 
comes gapless. The simplest picture of the pinball liquid 
is that the extra charge moves around like a " ball" avoid- 
ing the "pins" on A sublattice. The kinetic energy gain 
per single "ball" corresponds to the bandwidth of the 
honeycomb lattice which is obtained after depleting the 
A sublattice(pinned sites) from the original lattice. With 
increasing filling the energy gain will be affected by the 
correlation among "balls". We plot in Fig. (2J the esti- 
mation of Ei = (Eq — Ec)/N, where Eq is numerically 
obtained by the diagonalization of eq.QJ and the clas- 
sical binding energy is given by Ec = 3VN e , where N e 
denotes the electron number. The kinetic energy gain, 
1 23i I, increases with increasing number of filling, i.e. the 
number of "balls" that can move. In Fig.|2Ja) we fit the 
data by the power of 1/V as, E\ = c + C\/V + C2/V 2 , 
where Cj are fitting constants and cq gives the kinetic 
energy gain in the strong coupling limit. The obtained 
results at V — 00 are given as a function of filling in 
Fig- I3b) together with Ei's at finite V. In the strong 
coupling limit, E\(V — 00) seems to follow the band en- 
ergy of the honeycomb lattice at N e < 10 which supports 
the simplest pinball liquid picture we discussed above in 
the dilute case. As the filling approaches n = 1/2, E\ 
deviates from the honeycome one and saturates. This is 
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because the correlation among "balls" increases to avoid 
the loss of V as they becomes dense. Then the balls can 
no longer move independently with each other and the 
wave functions become relatively localized so that the 
kinetic energy gain do not increase much with N e . 

Although we see shortly that the pinball liquid is dom- 
inant over 1/3 < n < 1/2, its stability is still the 
most difficult to be guaranteed at half-filling because the 
largest number of contingents exists that might destroy 
the "pins" to some extent. We remind that the charges 
on B and C sites can hop by t to its neighboring site as 
in Fig. Of a) without the loss of binding energy when the 
neighboring population number of that charge remains 
unchanged. The example is given in Fig. Ofb) where 
some of the balls can move in the direction indicated 
by arrows. This rule also holds for the charges on A 
site in the irregular contingent states, e.g. when three 
neighboring hexagons have special kind of configuration 
as shown in Fig. Etc) . Fortunately, in this case the con- 
tingent cannot destroy the "pins" completely. This is 
because the charges which was originally a "pin" cannot 
go more than one site away without freezing the con- 
figuration of original two neighbors. Otherwise, one of 
these neighbors refills the original "pin" and the "balls" 
can move again smoothly. Therefore although the lo- 
cal fluctuation squeeze the amplitude of "pins" they re- 
main relatively stable and the collective destruction of 
"pins" does not take place. Considering the particle-hole 
symmetry, the same discussion holds by regarding the 
C-sublattice(hole) as "pins". Therefore, the "pins" and 
"balls" are not completely fixed and free, respectively, at 
half-filling. 

CORRELATION FUNCTIONS 

To further compare the stability of the pinball liquid 
under different fillings, we present the numerical estima- 
tion of order parameters for several number of charges 
N e =8-12 at N = 24, which correspond to the fillings, 
n=l/3-l/2. As discussed elsewhere, the pinball liquid is 
characterized by a three-body correlation 10]. Figure 01 
shows such function, Ps(i) = — n2)ni), at t = 1. 

As we see in Fig^Ja), the amplitude of "pins"(j4) near 
n=l/3 is almost 1 which decays only little with distance. 
This sturdily supports the strong coupling picture and is 
also consistent with the fact that the kinetic energy gain 
in Fig. Gla) at low fillings is almost identical with the 
honeycomb lattice band energy. With increasing N e the 
oscillation amplitude of P3 becomes suppressed and has 



larger decay. It is more clearly displayed in FigQJb) as 
a function of N e for several choices of V. However, the 
amplitude of "pins" are always larger for larger V, and 
we expect a long range order in the bulk limit at least in 
the very strong coupling region even at half-filling. 
SUMMARY 

We displayed the filling dependence of the new gap- 
less charge liquid denoted as pinball liquid, which was 
reported very recently at 1/3 < n < 1/2. Due to the 
particle-hole symmetry, the same results are obtained 
for 1/2 < n < 2/3 by just regarding holes as "pins". 
Despite the anticipation that the increasing number of 
particles towards half- filling might drive the "pins" to- 
wards destruction, the kinetic energy gain and the cor- 
relation function indicates that the pinball picture still 
holds while the balls become strongly correlated. 
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